Exploring Interatomic Coulombic Decay by Free Electron Lasers 
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To exploit the high intensity of laser radiation, we propose to select frequencies at which single- 
photon absorption is of too low energy and two or more photons are needed to produce states of an 
atom that can undergo interatomic Coulombic decay (ICD) with its neighbors. For Ne dimer it is ex- 
plicitly demonstrated that the proposed scheme to investigate interatomic processes by multiphoton 
absorption is much more efficient than with single-photon absorption of sufficiently large frequency 
as used until now. Extensive calculations on Ne dimer including all the involved nuclear dynamics 
and the losses by ionization of the participating states show how the low-energy ICD electrons and 
Ne + pairs are produced for different laser intensities and pulse durations. At higher intensities the 
production of Ne + pairs by successive ionization of the two atoms becomes competitive and the 
respective emitted electrons interfere coherently with the ICD electrons. It is also demonstrated 
that a measurement after a time delay can be used to determine the contribution of ICD even at 
high laser intensity. 

PACS numbers: 33.20.Xx, 32.80.Hd, 41.60.Cr, 82.50.Kx 



Inner- valence (IV) ionized states of isolated atoms and 
small molecules are usually electronically stable and may 
only relax via fluorescence decay or dissociation. When 
embedded in environment the IV ionized system can effi- 
ciently relax on a femtosecond time scale by transferring 
its excess energy to its neighbors and ionizing them |l| . In 
the last decade, many theoretical and experimental stud- 
ies of this process called Intermolecular or Interatomic 
Coulombic Decay (ICD) have been reported for various 
systems @, H| ■ The ultrafast relaxation driven by inter- 
atomic and intermolecular processes on the femtosecond 
time scale is of multidisciplinary importance in physics, 
chemistry, biochemistry, and biophysics [1 11 ■ 

There are many possibilities to initiate ICD in a sys- 
tem by single-photon absorption of sufficient energy to 
ionize an IV electron of the system. The new generation 
of light sources, free electron lasers (FELs), provide the 
opportunity to study radiation-matter interactions under 
extreme conditions, such as unprecedented high-intensity 
and ultrashort pulse durations. The presently operating 
FELs cover a wide range of short-wavelength radiation 
from VUV to soft X-ray (at FLASH facility [H) and 
up to hard X-ray (at LCLS facility [Hj]). Exposed to 
strong pulses, atomic and molecular systems may absorb 
a large number of photons creating differently charged 
ions and even bare nuclei 14{. Moreover, if the energy of 



a single-photon is insufficient to produce a selected state, 
the high-photon flux of a laser pulse enables the produc- 
tion of such states by multi-photon absorption [15!]. It is, 
therefore, rather natural to exploit these advantages of 
FELs for studying ICD processes. 

Very recently IH, 17 [, two possibilities to initiate ICD 
by multi-photon absorption in a cluster were discussed. 
In the collective ICD suggested in Ref. jjij two atoms 
are IV ionized by two photons in a cluster where a single 
IV vacancy cannot trigger ICD. The simultaneous relax- 



ation of the two IV vacancies can, however, undergo ICD 
by ionizing a third neighbor. In the process discussed in 
Ref. 17j, two or more atoms are excited and the relax- 



ation of one of the excited atoms results in the ioniza- 
tion of one of the other excited atoms. As discussed in 
Refs. 1(3, 17 1 , in intense fields these processes can provide 
efficient multi-photon ionization pathways in matter. 

In the mechanisms discussed in Refs. 3, 17 1 one can- 
not view the respective processes as the relaxation of 
an ionized or excited atom embedded in an environment 
as also a neighbor has to be ionized or excited. In the 
present work we intend to close the gap by investigating 
an atom which absorbs at least two photons and decays 
via ICD by ionizing its neighbor, providing a hitherto 
unexplored possibility to initiate ICD by an intense FEL 
pulse. Our results pave the way for further studies of in- 
teratomic processes by strong pulses of low-energy pho- 
tons. A single absorbed photon which ionizes an outer- 
valence (OV) electron of an atom (or molecule) embed- 
ded in environment does not trigger ICD. On the other 
hand, irradiating the system by an intense laser pulse 
with energy above the OV, but below the IV ionization 
threshold, can induce ICD. This particularly applies if 
the central frequency of the pulse is also tuned to the en- 
ergy difference between the IV and OV ionized states of 
the atom as these states will be strongly coupled by the 
intense field. For isolated atoms and diatomic molecules 
the coupling of two electronic states by an intense x-ray 
pulse of resonant frequency and the resulting Rabi oscil- 
lations of the population between these states have been 
recently studied 18-23j. 



The process proposed here can take place in any atom 
and molecule which possesses accessible states with suf- 
ficient excess energy above the ionization potential of its 
environment to allow for ICD to occur. However, we 
would like to suggest a candidate suitable for experi- 
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FIG. 1: (Color online) Ab initio PECs of the states relevant 
to the process {TJ of exposing the Ne dimer to an intense laser 
pulse. Lower panel: Of the Ne2 ground state. Middle panel: 
PECs of the Ne + (2p _1 )Ne ionic OV states. Upper panel: 
PECs of the Ne+(2s _1 )Ne ionic IV states which decay by ICD 
and of the dications Ne + (2p _1 )-Ne + (2p _1 ) which are the fi- 
nal states of ICD. The photon frequency of u = 26.888 eV 
used in the calculations is indicated in the figure. Note the 
different energy scales in all panels. 



mental verification and at the same time this candidate 
should be amenable to high quality calculations includ- 
ing the nuclear dynamics of the whole process. The Ne 
dimer fulfills these prerequisites. We report computed 
data and ICD electron spectra which may directly be 
compared with experiments. Such an experiment to ver- 
ify the proposed mechanism is feasible at present, i.e. at 
the FLASH facility. 

For the example chosen the process can be schemati- 
cally described as follows: 
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The field ionizes a 2p electron of a Ne atom in the dimer 
and couples resonantly the resulting OV ionic states with 
the respective IV ionic states. The vertical double-arrow 
in (JTJ indicates the strong coupling between the OV and 
IV ionic states. Fig. [1] depicts the ab initio computed 
potential energy curves (PECs) of the states relevant for 
the process. At the equilibrium internuclear distance 
(R e = 3.1 A) of the Ne dimer, the OV ionic states 
Ne + (2p _1 )Ne (middle panel) are nearly degenerate and 
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FIG. 2: (Color online) Integral yield of the Ne + ion pairs 
produced by exposing Ne dimers to a Gaussian laser pulse of 
60 fs duration as a function of the peak intensity of the pulse. 
Open circles: Yield computed for the process {TJ for a pulse 
with a central frequency uj = 26.888 eV. At least two photons 
are required for the process. Open squares: Yield computed 
for a pulse with a central frequency of of ui = 53.776 eV. The 
ICD is triggered directly by the absorption of one photon. At 
the intensities shown the ICD initiated by process ([T]) via two- 
photon absorption is clearly the more efficient mechanism for 
the production of these ionic fragments. 



this holds also for the IV ionic states Ne + (2s _1 )Ne (up- 
permost panel). The near degeneracy becomes rather 
perfect at larger internuclear distances. To resonantly 
couple the OV and IV states by an intense laser field we 
choose a pulse with a central frequency of u) — 26.888 eV, 
see Fig.[TJ which is the energy difference between the sta- 
tistically weighted average of the two groups of computed 
energies at 3.1 A. The IV ionic states shown in the up- 
per panel of Fig. [T] decay by ICD producing the strongly 
repelling Ne + (2p _1 )-Ne + (2p _1 ) states (also shown in 
this panel) which undergo a Coulomb explosion and low- 
energy ICD electrons. Importantly for experiment, the 
photoelectrons produced in the first step of {1} initiating 
the whole process have an energy of around 5.32 eV which 
is well separated from the energy of the ICD electrons. 
Finally we stress that Ne+ (2p~ 1 )-Ne+ (2p~ 1 ) can be com- 
petitively produced by the direct ionization of the neu- 
tral Ne atom in Ne + (2p _1 )Ne produced in the first step 
of fl} and similarly Ne + (2s _1 )-Ne + (2p _1 ) (not shown in 
Fig. [I} by the direct ionization of Ne + (2s _1 )Ne resulting 
in the second step of ([T]). All of these states and processes 
are included in the computations. 

In order to compute the process (JTJ) we combine the 
previously developed theoretical and computational ap- 
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proaches 

effect in intense laser fields 



to evaluate the resonant Auger decay 
(see [H for details). The 



nuclear dynamics have been calculated employing the 
efficient Multi-Configuration Time-Dependent Hartree 
(MCTDH) method and code [HEl. Thereby we have 
utilized the ab initio PECs and the ab initio ICD tran- 
sition rates reported in Ref. 27] . The values of the 



electron transition matrix elements were extracted from 
the experimental photoionization cross section of the 
Ne atom (a 2p = 7.8 Mb at 28.4 eV [11]) and the ex- 
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FIG. 3: (Color online) Electron spectra after exposure of Ne2 
to a 60 fs Gaussian pulse for three intensities (solid curves, 
'Total'). For comparison the individual contributions of the 
ICD process (broken curves, 'ICD') are also shown. For each 
intensity two numbers are reported in brackets: the fraction of 
neutral dimers which have survived the pulse, and the fraction 
of dimers which have become singly ionized. 

perimental 2s _1 —> 2p -1 radiative decay rate of Ne + 
(r r = 4.8 fieV [29[ corresponding to the radiative life- 
time of r r ~ 0.14 ns). 

Fig. [5] depicts the integral yield of the Ne + ion pairs 
produced by exposing the Ne dimer to a Gaussian-shaped 
pulse of 60 fs duration as a function of the pulse peak 
intensity. For comparison, we remind that the ICD life- 
time is ~80 fs [27| ■ The yield computed for the presently 
studied process ([T]) is depicted by open circles. At the rel- 
atively low intensities considered in the figure, the ICD 
initiated by two-photon absorption is by far the dominant 
mechanism for the production of these ionic fragments. 
As expected, the yield grows quadratically with the field 
intensity as the fit in Fig. [2] to the five points lowest in 
intensity shows. How does this yield compare with the 
yield obtained in the traditional way of initiating ICD 
by the absorption of a single-photon of energy above the 
IV ionization threshold? To answer this basic question 
we also computed the ICD triggered directly by a strong 
pulse with twice the central frequency. The results are 
shown as open squares in the figure. The computed yield 
grows linearly with the intensity, as the fit in Fig. [2] to 
the first five points lowest in intensity indicates. Inter- 
estingly, the process (JTJ is much more efficient at the 
considered intensities, although it requires the absorp- 
tion of at least two photons. This finding gives further 
motivation to investigate interatomic and intermolecular 
processes by intense lasers. 

At intensities larger than 10 11 W/cm 2 the computed 
integral yield starts to deviate from the fit. Indeed, as 
the field intensity grows other mechanisms than ICD in- 
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FIG. 4: (Color online) Electron spectrum after exposure of 
Ne2 to a 30 fs Gaussian pulse of intensity 1 x 10 13 W/cm 2 
(solid curve, 'Total'). Also shown are the individual contri- 
bution of the ICD process (dashed curve, 'ICD'), the spec- 
trum measured 30 fs after the pulse maximum (dotted curve, 
'Totabjo'), and the difference between the 'Total' and the 
'Totabjo' spectra (dash-dotted curve). 



duced by the field start to play an important role. As 
demonstrated in Refs. [19l421j. the (direct) ionization of 
all the states participating in the ICD process results 
in losses (leakages) of population of these states. Par- 
ticularly relevant and interesting is the direct ionization 
of Ne + (2p _1 )Ne mentioned above. This ionization and 
the ICD populate the Ne + (2p" 1 )-Nc+(2p" 1 ) states co- 
herently, thus naturally inducing interference effects in 
the resulting electron spectra 12 1 |. As has been dis- 
cussed and demonstrated in Refs. |22ll23j|. an intense field 
also gives rise to light-induced non-adiabatic effects. All 
these mechanisms were simultaneously incorporated in 
the present calculations. 

Let us now turn to the spectra of the electrons emit- 
ted in the process (|T|). The electron spectra computed 
for Gaussian-shaped 60 fs pulses of three different peak 
intensities are depicted in Fig. [3] by solid curves (la- 
beled 'Total'). Although the intensities chosen are not 
the strongest FEL pulse intensities available at present, 
they cannot be accessed by any conventional synchrotron 
sources. After the pulse is over and all the IV states 
populated by it have decayed via ICD, a fraction of the 
neutral dimers, of course, survives and remains in their 
ground state. For each intensity this fraction is reported 
in the figure (first number in the bracket below the indi- 
cated intensity). Also of interest are the corresponding 
fractions of the dimers which have become singly ionized 
NeJ dimers. These numbers are reported in the figure 
next to the neutral fraction. At the lowest considered 
peak intensity (5 x 10 11 W/cm 2 ), the pulse manages to 
ionize 13% of the neutral dimers. Of these, more than 
half (7.2%) remain as NeJ ions and the rest (5.8%) par- 
ticipate further in the ICD. Even at highest considered 
intensity of 5 x 10 12 W/cm 2 , the ionization of the ground 
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FIG. 5: (Color online) Electron spectrum after exposure of 
Ne2 to a 30 fs Gaussian pulse of intensity 5 x 10 12 W/cm 2 
(solid curve, 'Total'). For analysis the individual contribu- 
tion of the ICD process (dashed curve, 'ICD'), the individ- 
ual contribution of the direct ionization, i.e., excluding ICD 
from the calculations (dotted curve, 'Direct'), and the results 
obtained excluding only the ionization of Ne + (2s _1 )Ne to 
produce Ne + (2s _1 )-Ne + (2p _1 ) (dashed-dotted curve, 'Sub- 
total') are shown. 



and of the ionic states are still far from being saturated: 
about a quarter of the neutral dimers survive, and of the 
three quarters which have been ionized about 30% remain 
as NeJ ions and 45% have undergone further ionization 
and have emitted an electron which contributes to the 
spectrum shown. The calculations indicate that sequen- 
tial single-photon ionization dominates over simultaneous 
two-photon absorption. 

As discussed above, not all of the low-energy electrons 
in the computed spectra in Fig. [3] are produced via the 
ICD. In order to reveal the individual contribution of the 
ICD process, we performed additional calculations. In 
these all the pathways to produce Ne + -Ne + dicationic 
states by direct ionization of Ne^~ ions were excluded, 
but all losses from the neutral and ionic states were still 
taken into account. The results of these approximate 
calculations are depicted in Fig. [3J by broken curves (la- 
beled TCD'). The difference between the solid and broken 
curves at each peak intensity represents the contribution 
to the spectrum deriving from the direct ionization path- 
ways and the interference between the direct and ICD 
pathways. 

As seen in Fig. [3J almost all low-energy electrons 
are produced via the ICD for intensities below 5 x 
10 11 W/cm 2 . Interestingly, the shape of the ICD spec- 
trum computed for this peak intensity is rather similar 
to that of the experimental ICD spectrum measured in 
Ref. Q by synchrotron radiation in the single-photon ab- 
sorption mode. ICD is still by far the dominant pathway 
for production of low-energy electrons at 1 x 10 12 W/cm 2 
(middle panel), but the direct ionization pathways be- 
come more and more relevant as the intensity increases. 



This leads to a change of the shape of the electron spec- 
trum, indicating already now (see below) that the direct 
pathways preferentially produce electrons at somewhat 
higher energies than ICD. 

How to identify at strong field the individual contribu- 
tion of ICD to the spectrum? We may exploit the fact 
that the two kinds of pathways for the production of the 
final ionic fragments take place on different timescales. 
Obviously, the direct ionizations are operative only when 
the pulse is on. The ICD process itself is, in turn, inde- 
pendent of the pulse duration and in the present example 
the ICD lifetime is ~ 80 fs 27]. In order to separate the 



two kinds of pathways, one has to apply shorter pulses 
and invoke measurements in the time domain. 

To demonstrate this idea we show in Fig. 2] electron 
spectra computed for a Gaussian pulse of half the dura- 
tion (30 fs) and twice the intensity (1 x 10 13 W/cm 2 ) of 
the pulse used to obtain the spectrum in the upper panel 
of Fig. [3j As usual, the total spectrum (solid curve) is 
obtained at long time after the process is over. Measur- 
ing the spectrum after a time delay of 30 fs from the 
pulse maximum, i.e., just after the pulse expired, pro- 
duces the spectrum labeled 'Totai3o' depicted as a dotted 
curve. During the time the pulse is on the direct ioniza- 
tion is expected to strongly dominate the production of 
the spectrum and the difference between the 'Total' and 
the 'Totalso' spectra is thus expected to reflect the indi- 
vidual contribution of the ICD. This difference, shown in 
Fig.[4]as a dash-dotted curve, is indeed very similar to the 
spectrum computed excluding the direct ionization path- 
ways (dashed curve). This finding confirms that most 
ICD electrons are emitted after the pulse has expired. 

It is evident from Figs. [3] and 0] that direct photoioniza- 
tion and ICD preferentially produce electrons in different 
energy ranges. This finding is further analyzed in Fig. [5] 
for a Gaussian pulse of 30 fs duration and 5 x 10 12 W/cm 2 
peak intensity. The individual contribution of the di- 
rect ionization is shown by the dotted curve (labeled 'Di- 
rect'), and that of ICD by the dashed curve. Clearly, ICD 
mainly contributes to the electron energy range 0-1 eV, 
and the direct channel produces electrons mainly between 
0.5 and 1.5 eV. Most importantly, the two contributions 
are by no means additive and the strong interference be- 
tween the pathways determines the final shape of the 
spectrum. As a last point we mention that at high inten- 
sities also the ionization of Ne + (2s _1 )Ne which can decay 
by ICD contributes and cannot be neglected. At least a 
third photon must be absorbed in the process for this 
ionization which produces Ne + (2s _1 )-Ne + (2p -1 ). The 
impact of this ionization can be seen in Fig. [5] (dash- 
dotted curve labelled 'Sub-total'). 

In conclusion, a general mechanism to initiate ICD by 
an intense laser pulse is proposed. The initiation requires 
the absorption of at least two-photons. Nevertheless, it 
is found here to initiate interatomic processes much more 
efficiently than by utilizing traditional one-photon ioniza- 
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tion schemes with photons of sufficiently high energies. 
The process is exemplified by extensive calculations on 
the Ne dimer which take account of losses via ionization 
induced by the intense field and incorporate the quantum 
nuclear dynamics on all the participating electronic levels 
including their coupling by the pulse. At moderate field 
intensities, ICD is the main mechanism for the produc- 
tion of low-energy electrons. At higher intensities other 
ionization mechanisms become relevant as well. Impor- 
tantly, they coherently interfere with ICD enriching the 
process. Since these other mechanisms are only operative 
while the pulse is on, there is an interesting interplay be- 
tween them and ICD which can be controlled by varying 
the pulse. In particular, pump-probe techniques can be 
used to identify the contribution of ICD to the electron 
spectrum even at higher intensity by a single additional 
measurement after a time delay. It is our opinion that the 
conclusions drawn are rather general and apply to many 
systems. We hope that the present results will stimulate 
experiments. 

The research leading to these results has received fund- 
ing from the ERC under the EU's FP7, AIG No. 227597. 
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I. THEORY: NUCLEAR DYNAMICS HAMILTONIAN 



In this supplementary information we provide an extended representation of the presently 
applied theoretical approach. It implies a general formulation of the time- dependent theory 
for the nuclear wave packet propagation (see, e.g., Refs. [S1-S5]), which has recently been 
extended and applied to evaluate the resonant Auger decay effect of atoms and molecules 
in intense laser fields [S6-S8] . The extended approach and all necessary derivations of the 
theory can be found in Refs. [S5-S8]. In order to make it easier for the interested reader 
to follow our approach without restoring to the literature, its essential points and 'working' 
equations describing the presently studied process are collected and discussed below. 

We start with the schematic representation of the example chosen in the main text: 

Ne 2 Ne + (2j 5 - 1 )Ne + e ph 

Ne + (2s- 1 )Ne (SI) 

J, ICD 

Ne+(2p- 1 ) + Ne+(2p- 1 ) + e ICD 

The strong laser field ionizes a 2p electron of a Ne atom in the dimer populating outer- valence 
(OV) ionized states Ne + (2p _1 )Ne and producing a photoelectron e ph . The same pulse with 
the resonant central frequency u = 26.888 eV couples also the resulting OV and inner-valence 
(IV) ionized states Ne + (2s _1 )Ne, as indicated by the vertical double arrow in Eq. (jSip . The 
IV ionic states of the dimer relax by interatomic Coulombic decay (ICD) producing the 
strongly repelling Ne + (2p _1 ) — Ne + (2p~ 1 ) states which undergo a Coulomb explosion and 
low-energy ICD electrons eico- The Ne + (2p _1 ) — Ne + (2p _1 ) states can be competitively 
produced by the direct ionization of the neutral Ne atom in Ne + (2p _1 )Ne states produced 
in the first step of fISip . Similarly, the ionization of IV ionized states Ne + (2s _1 )Ne resulting 
in the second step of ( 1ST]) to directly produce the Ne + (2s _1 ) — Ne + (2p _1 ) states cannot be 
neglected at high intensities of the field. All of these states and processes are included in 
the present computations. 

In order to theoretically describe the dynamics of the process ( ISljl as a function of time, 
we solve the time-dependent Schrodinger equation for the Ne 2 and its interaction with the 
field (atomic units e = m e = h = 1 are used throughout) 

i4?(t) = #(*)*(*) = {H nuc + Hei + DE{t)) tt(t). (S2) 

where D represents the dipole transition operator. We assume a coherent and monochro- 
matic Gaussian-shaped pulse of duration r centered at to {git) = e - ^ - * ^ 2 / 1 " 2 ) with linear 
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polarization of the field along the z axis and central frequency u: 



£{t) = £o(t) cosut = £og{i) cosut. (S3) 

Here £q is the peak amplitude, and the pulse-shape function g(t) varies slowly on the 
timescale of 2rr/u. The cycle-averaged intensity of the field is given in atomic units via 
(1 a.u. = 6.43641xl0 15 W/cm 2 ) 

I(t) = -^-{£ g(t)} 2 , (S4) 

where a = 1/137.036 is the fine structure constant. 

Following Refs. [S5-S8] the total wave function of a rotating molecule as a function of 
time can be represented via the ansatz including all electronic states participating in the 
process (jSip . It includes: the ground state of the Ne2 with the electronic wave function $/; all 
intermediate Ne + (2p _1 )Ne (OV-ionized) and Ne + (2s _1 )Ne (IV-ionized) states plus outgoing 
photoelectron of energy e p h with the total electronic wave functions and ; and all 
final Ne + (2p _1 ) — Ne + (2p _1 ) (OVOV doubly-ionized) states plus outgoing photoelectron and 
ICD electron of energy Ejcd with the total electronic wave function ^ S Qy^f 

*(f) = MWi + E / ^ov i (e ph ,t)^ J de ph + Y / J ^i Vj (e ph ,t)^.de ph + 

3 j 

^ovoVj( £ ph, £icd, t)^ovov° d £ phdejcD, (S5) 

where index j numerates all possible states within each subset of states. The functions 
\&j(t), ^oVj(s P h,t), $ 'iVj(£ph,t) and ^ovov 3 { £ P h, £ icd, t) are the time-dependent wave pack- 
ets propagating on the potential energy surfaces of the included electronic states. It should 
be remembered that these wave packets depend on the nuclear coordinates R and 9. For 
brevity, we will explicitly show the dependence on the nuclear coordinates only in the final 
Hamiltonian matrix (lS8p governing the nuclear dynamics of the process (ISip . 

In order to obtain the final set of equations for the propagation of the nuclear wave 
packets, we substitute the ansatz flS5|) in the time-dependent Schrodinger equation ({32]) and 
project the result onto each electronic state. As explicitly demonstrated in Refs. [S5-S8], 
in the derivation we imply the rotating wave approximation [S9], the local approximation 
[S10,S11], and redefine ('dress') the time dependent wave packets as follows 

*o Vj {e ph , t) = V OVj {e ph , t) e +iut (S6a) 

(e P h, t) = %i Vj (eph, t) e +2lult . (S6b) 
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^OVOVj ( £ ph, £lCD, t) — ^OVOVj ( e ph, &ICD, t) e +2lult (S6c) 

Of course, all electronic states within the subsets of IV and OV singly-ionized, and OVOV 
doubly-ionized final states shown in Fig. 1 of the main text participate in the dynamics of the 
considered process and have to be included in the ansatz (1S5I) . For the sake of transparency 
of presentation we explicitly write in the ansatz (1S5I) below only one electronic state from 
each of the subsets. In particular, we show only the Ne + (2p _1 )Ne OV-ionized state of 2 U g 
symmetry coupled by the field with the Ne + (2s _1 )Ne IV-ionized state of 2 £+ symmetry, 
which subsequently decays via ICD into one of the OVOV doubly-ionized final states. The 
extension of the expressions to the case of all participating states is straightforward. 

Let us collect the individual nuclear wave packets contributing to the total wave function 
(155]) into a single vector 

/!*/(*)> 

|*ov(£pfc,t)) 
\Viv(£ph,t)) 

t))J 



\y(e ph ,e ICD ,t)) 



(S7) 



and define the matrix Hamiltonian 



K(R,9,t) = T(R,6)+ 



( Vj(R) - |rf (t) 







\ 



d x (t) sin^+ 
+d z {t) cos 9 



Vov(R) - |rgvw- 

+e ph - u 



(Dl(t)-iW*{t))Bm6 







(D x (t)-iW(t))sm8 



V IV {R)+e ph -2u- 







d x (t) sin 6+ 
+d z {t) cos 9 



V ICD (R) 



Vovov(R)+ 
+e ph + e ICD -2u / 



(S8) 



With these notations, the final set of equations describing the propagation of the nuclear 
wave packets (IS 7ft now takes on the following compact form 



\^(e ph ,e IC D,t)) = H(t) \^{e ph ,e IC D,t)). 



(S9) 



The matrix H can be viewed as the effective Hamiltonian governing the two-dimensional 
nuclear dynamics in the process (IS II) in an intense laser field. Below we summarize the 
physical meaning of each term of the Hamiltonian matrix flSSj) and provide their explicit 
expressions. 
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T(R, 9) is the common nuclear kinetic energy operator for the vibrational motion along 
the internuclear distance R and the rotational motion described by the angle 9 between the 
polarization vector of the laser pulse and the molecular axis. The functions Vi(R), Vov(R), 
Viv(R), and Vovov(R) on the diagonal are the potential energies of the corresponding 
electronic states. These potential curves of the Ne2 were computed ab-initio in the present 
work as described in Refs. [S12,S13] and are depicted in Fig. 1 of the main text. The 
couplings between the various electronic states are given by the following matrix elements 
of the total Hamiltonian (IS2I) 

($ovl-fr(*)l$J> = [d x sm6 + d z cos9] £°|21 e -^ = [d x (t) sin 9 + d z (t) cos 9} e~ iw \ (SlOa) 
(^% h \H(t)\^) = D x sm9^^-e-^ = D x (t) sin9 e*> \ (SlOb) 

($ovov D l#(*)l$o£> = sin d + d ~* cos 9 e~ lujt = d x {t) sin 9 + d z {t) cos 9 e _<w *, 

(SlOc) 

^Tvov D mm7v) = v icd (r), (siod) 

In Eqs. (1S10al - [Sl0cl) . the rotating wave approximation [S9] has already been utilized, and, 
in contrast to the rapidly oscillating factor e~ luJt , the functions D(t) and d(t) vary slowly on 
the timescale of 2ir/u. Let us discuss these transition matrix elements. 

The matrix element (ISlOaj) describes the ionization of a 2p electron of a Ne atom of the 
dimer in the ground electronic state via the x and z components of the dipole transition 
operator to produce OV-ionized state Ne + (2p _1 )Ne and a photoelectron e p h- The matrix 
element fISlObj) represents the excitation of a 2s electron of Ne + ion into the vacant 2p orbital 
via the x component of the dipole transition operator and couples, thereby, the selected 
OV 2 II 9 and IV 2 E+ states of the singly-ionized dimer. The matrix element (ISlOcj) is similar 
to (ISlOaj) and corresponds to the ionization of the neutral Ne atom in the Ne + (2p _1 )Ne state 
to directly produce the Ne + (2p _1 ) — Ne + (2p _1 ) final doubly-ionized state and ICD electron 
Eicd- Eq. fISlOdj) defines the matrix element for the ICD transition, i.e. of the decay of the 
selected IV-ionized state into the OVOV doubly-ionized final state with emission of the ICD 
electron via the Coulomb operator (note its explicit dependence on the internuclear distance 
R [S14,S15]). 

The ionization of the neutral dimers in the first step of the process (IS 1 j) transfers the 
population from the ground electronic state to the OV-ionized states of the dimer by the 
transition matrix element (ISlOaj) (matrix element if 21 in (1S8j) ). The leakage of the corre- 
sponding population from the ground state due to direct photoionization into all possible 
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final ionic states is described by the imaginary time- dependent term — ^T ph (t) on the re- 
spective diagonal (Hu) in the Hamiltonian ( 1S8|) . Its explicit expression was obtained in [S6] 
in the local approximation and reads: 

rf(t) = 2vr^|< /z (t)| 2 , (Sll) 

3 

where d 3 x , x (t) are the dipole transition matrix elements (ISlOaj) for the direct ionization of 
the ground electronic state into all possible final ionic states numerated by superscript j. 
The total probability for the direct photoionization of the ground state (IS 1 1 H is identical to 
the quantity 7 ph (t) introduced in Ref. [S16]: 

rf (t) = Y h (t) = of i{t)/u, (Si2) 

where a ph is the total direct photoionization cross section of the ground state for the exciting- 
photon energy oj, I(t) is the field intensity (154)) . and the quantity I(t)/u stands in Eq. (IS 1 2|) 
for the photon flux. 

Owing to the redefinition (1561) . the potential energies of the IV and OV ionized states 
with a photoelectron on the respective diagonal elements of H are 'dressed' by the field (i.e., 
Vov(R) + £ph ~ w an d Vjy(R) + £ P h — 2w). They are also augmented by time-dependent 
imaginary terms — ^T P Q V (t) and — jT p IV (t), respectively, similar to Eq. (1S12|) . which describe 
the leakages of the populations of the OV and IV ionized states due to direct ionizations of 
the neutral Ne atom in the Ne + (2£ _1 )Ne to produce all possible final Ne + (2£ _1 ) — ~Ne + (2p~ 1 ) 
dicationic states. Both T P Q V {t) and T P y{t) probabilities can be computed by equations like 
Eq. (lS12j) . The time- independent imaginary term — ^T I I y D {R) represents the losses of the 
population of the initial ICD state (i.e., IV-ionized state Ne + (2s _1 )Ne) by ICD transition, 
i.e., the total ICD rate. This rate is given in the local approximation by the Coulomb matrix 
element f jSlOdj) as 

r^( J R) = 2 7 r^|\// CD (i?)| 2 , (S13) 



where the summation over index j runs over all possible ICD channels 

As was demonstrated in Ref. [S6], the usual direct coupling D x {t) between the IV and 
OV ionized states through the laser field (by the 2s 2p excitation-deexcitation of Ne + in 
the present case) is augmented by an additional time-dependent term — ^W(t). This term 
appears only if the photoionization from the OV-ionized state and ICD transition from IV 
ionized state are simultaneously treated as required. Explicitly, it reads [S6] 

W(t) = 2nY,dl(t){Vi CD y , (S14) 
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where the summation over j runs over all possible final dicationic states accessible by both 
channels. The whole coupling between the 'dressed' OV and IV ionized states is non- 
Hermitian and only operative as long as the pulse is on. The forth row of the Hamiltonian 
matrix (158]) describes the nuclear dynamics on the Ne + (2p _1 ) — Ne + (2p _1 ) final OVOV- 
ionized state produced by the emission of the ICD electron Eicd- It is populated coherently 
by the ICD transition from the IV- ionized state (matrix element (jSlOdj) ) and by direct 
photoionization from the OV-ionized state (matrix element (ISlOcj) ) . all at a given kinetic 
energy Sjcd of the emitted electron. The created wave packet propagates on the 'dressed' 
potential energy surface Vovov(R) + e ph + £/cd — 2a;. 

The nuclear wave packet ^ovovi^ph, £icd, t) contains the information on the coincident 
photoelectron and ICD electron spectrum. The coincident spectrum can be computed as 
the norm of ^ovov(^ph,£icD,t) at long times [S5]. Integration of the coincident spectrum 
over the energy of the photoelectron yields the individual ICD spectrum computed in the 
present study via 

v{sicd) = lim / (^ovoVji^ph^iCD.t^ovov^ph^iCD.t)) de ph . (S15) 

As a last point we mention that at high intensities also the ionization of the initial ICD states 
Ne + (2s _1 )Ne to directly produce the Ne + (2s _1 ) — Ne + (2p _1 ) states contributes to the final 
spectrum. To this end at least a third photon must be absorbed in the process flSlj) . In order 
to incorporate this mechanism in the calculations one has also to include the corresponding 
electronic states in the ansatz (1S5[) which is straightforwardly done. 

As has been demonstrated in Refs. [S7,S8], the 'dressed' states of a diatomic molecule can 
exhibit intersections of the corresponding two-dimensional complex potential energy surfaces 
in the R and 9 space. The dynamical variables R and 9 enter the Hamiltonian matrix (]S8j) 
explicitly. In the present case of ICD, the field couples the 'dressed' singly- ionized electronic 
states of the Ne^ (i.e., IV and OV states as discussed above). The rotational degree of 
freedom is involved in the nuclear dynamics only due to the presence of the laser field (the 
coupling matrix elements H32 and -£^3 are proportional to sin 9) and only when the pulse 
is on. Very important, the non-adiabatic couplings, i.e., matrix elements of the nuclear 
momenta along R and 9, between the two 'dressed' electronic states are singular at these 
intersections [S17]. This gives rise to dramatic dynamical effects [S7,S8]. 

Without the imaginary terms in the Hamiltonian C IS8 j) . the two-dimensional potential en- 
ergy surfaces Vov{R) and Viv(R) — u) of the 'dressed' electronic states exhibit an intersection 
at 9 = (where sin# = 0), which is a conical intersection [S18-S20]. Due to the presence of 
the ICD width and leakages by ionization in (158]) . the situation becomes more complicated. 
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The two potential energy surfaces in R and 9 space obtained by diagonalizing the electronic 
Hamiltonian H(i2, 9,t) — T(R,6) in Eq. (1S8|) are now complex and generally exhibit two 
intersecting points at which the real as well as the imaginary parts of the two electronic 
energies become degenerate [S21]. This analogue of a conical intersection in the continuum 
has been named doubly intersecting complex energy surfaces (DICES). The impact of the 
light-induced DICES on the nuclear dynamics accompanying resonant Auger decay of the 
core-excited CO molecule in intense laser fields has recently been studied in Ref. [S8]. The 



non-adiabatic effects of DICES are naturally incorporated in the Hamiltonian matrix (JS8J) 
and, thus, were included in the present study of ICD. 

To be able to carry out the two-dimensional calculations on the coupled complex energy 
surfaces we employed the efficient Multi-Configuration Time-Dependent Hartree (MCTDH) 
method [S22] and code [S23]. The theoretical study of the nuclear dynamics problem requires 
the potential energies of the electronic states participating in the process and electronic 
transition rates between. In the calculations we have utilized the presently computed ab 
initio energy curves (see Fig. 1 of the main text) and the ab initio ICD transition rates 
reported in Ref. [S24] . The values of the electron transition matrix elements were extracted 
from the experimental photoionization cross section of the Ne atom (o"2 P = 7.8 Mb at 28.4 eV 
[S25]) and the experimental 2s _1 — > 2p~ 1 radiative decay rate of Ne + (r r = 4.8 /ieV [S26] 
corresponding to the radiative lifetime of r r ~ 0.14 ns). 
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